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Abstract. It is found that the finite-energy spectral properties of the one- 
dimensional Hubbard model are controlled by the scattering of charged 77-spin-zero 
2z/-holon composite objects, spin-zero 2^-spinon composite objects, and charged rj- 
spin-less and spin-less objects, rather than by the scattering of independent 77-spin 
1/2 holons and spin 1/2 spinons. Here v = 1,2,.... The corresponding S matrix is 
calculated and its relation to the spectral properties is clarified. 
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The description of the microscopic scattering mechanisms behind the unusual finite- 
energy spectral properties observed in low-dimensional materials remains until now 
an interesting open problem. The one- dimensional (ID) Hubbard Hamiltonian is the 
simplest model for the description of electronic correlations in a chain of N a sites. It 
reads H = f + UD-[U/2][N-N a /2], where f = -iE ff=u Ej"i[ c L c J+V + tc -] is the 
kinetic- energy operator, D = J2j n j, \ the electron double-occupation operator, N = 
J2j.a^j,(T the electron number operator, and the operator ct creates a spin-a electron 
at site j. In contrast to other interacting models pQ and in spite of the model exact 
solution [2], until recently little was known about its finite-energy spectral properties 
for finite values of the on-site repulsion U. Recently, the problem was studied by the 
pseudofermion dynamical theory (PDT) introduced in Refs. 011], whose predictions 
agree quantitatively for the whole momentum and energy bandwidth with the peak 
dispersions observed for the TCNQ stacks by angle-resolved photoelectron spectroscopy 
in the quasi- ID conductor TTF-TCNQ and are consistent with the phase diagram 
observed for the (TMTTF^X and (TMTSF^X series of compounds More recently, 
results for the TTF-TCNQ spectrum consistent with those of the PDT were obtained 
by the dynamical density matrix renormalization group method [6,. Within the PDT, 
the finite-energy spectral properties are controlled by the functional character of the 
pseudofermion anticommutators jJE]. However, the relation of these anticommutators 
to the elementary-excitation S matrix remains an open question. Moreover, the fact 
that these anticommutators do not couple quantum objects with different r^-spin or 
spin projections seems to be inconsistent with the form of the S matrix for elementary 
excitations calculated in Refs. [3 IE]- Thus, the study of the relation of the PDT to the 
elementary-excitation scattering is an important issue both for the clarification of that 
apparent inconsistency and the further understanding of the scattering mechanisms that 
control the exotic finite-energy spectral properties of low-dimensional materials and of 
the new quantum systems described by cold fermionic atoms on an optical lattice Ell- 
in this Letter the above problems are solved by identifying the active scatterers and 
scattering centers which control the dynamical properties, calculating their S matrix, 
and clarifying its relation to the spectral properties. Moreover, the connection to the 
S matrix of Refs. [3 Ej is also clarified. The number of lattice sites N a is considered 
large, units of Planck constant and lattice spacing one are used, and the lattice length 
is denoted by L = N a and the electronic charge by — e. The densities n = N/L and 
spin densities m = [N^ — N^\/L are in the domains < n < 1 and < m < n, 
respectively. The above Hamiltonian commutes with the generators of the r^-spin and 
spin SU{2) algebras [TH]. Here the ?7-spin and spin values of an energy eigenstate are 
called 7] and S, respectively, and the corresponding projections r\ z and S z . A key result 
needed for our study is that all energy eigenstates of the model can be described in 
terms of occupancy configurations of r^-spin 1/2 holons, spin 1/2, spinons, and ^-spin- 
less and spin-less cO pseudoparticles [TU]. Below, the notation ±1/2 holons and ±1/2 
spinons is used according to the values of 77-spin and spin projections, respectively. The 
electron - rotated-electron unitary transformation JU] maps the electrons onto rotated 
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electrons such that rotated-electron double occupation, no occupation, and spin-up and 
spin-down single occupation are good quantum numbers for all values of U. The ±1/2 
holons of charge ±2e and zero spin and the charge-less ±1/2 spinons are generated 
from the electrons by that unitary transformation. The corresponding holon and spinon 
number operators M c , ±1/2 and M S) ±i/2, respectively, are of the form given in Eq. (24) 
of JU] and involve the electron - rotated-electron unitary operator. While the —1/2 
and ±1/2 holons refer to the rotated-electron doubly occupied and unoccupied sites, 
respectively, the —1/2 and ±1/2 spinons correspond to the spin degrees of freedom 
of the spin-down and spin-up rotated-electron singly occupied sites, respectively. The 
charge degrees of freedom of the latter sites are described by the spin-less and ^-spin- 
less cO pseudoparticles, which are composite objects of a charge — e chargeon and a 
charge ±e antichargeon ^OJ- The cv pseudoparticles (and sv pseudoparticles) such that 
v = 1,2,... are 77-spin singlet (and spin singlet) 2z/-holon (and 2z/-spinon) composite 
objects. Thus, M a> ±i/ 2 = £ a ,±i/2±IZ£Li v N au where a = c, s, N au denotes the number 
of av pseudoparticles, and L C) ±u 2 = r\ ± r\ z and L Sj ±i/ 2 — S ± S z gives the number of 
±1/2 Yang holons and ±1/2 HL spinons, respectively. Those are the holons and spinons 
that are not part of composite pseudoparticles. All energy eigenstates can be described 
by occupancy configurations of av pseudoparticles, —1/2 Yang holons, and —1/2 HL 
spinons [TU]. For the ground state, N c0 = N, N s i = iVj, N au = L Cj _x/ 2 = L S) _x/ 2 = 
for av 7^ cO, si. 

In our study we consider the pseudofermion subspace (PS), which is spanned by the 
initial ground state \GS) and all excited energy eigenstates contained in 0\GS), where 
O is any one-electron or two-electron operator. In reference [3] it is shown that within 
the PS there is a unitary transformation that maps the av pseudoparticle or hole onto 
the av pseudofermion or hole, respectively. These objects differ only in the discrete 
momentum values. The av pseudoparticle or hole has discrete bare-momentum values 
qj = [2ii/L}Ij' u such that are consecutive integers or half-odd integers jU3j. These 
values are good quantum numbers whose allowed occupancies are one (pseudoparticle) 
and zero (hole) only. The av pseudofermion or hole has discrete canonical-momentum 
values given by, 

Qj = Q(0j) = Qj + Qtu(Qj)/ L , (!) 

where j = 1,2,...,N* V , N* u = N au ± N^ v , and N^ v denotes the number oi av 
pseudofermion holes, which equals that of av pseudoparticle holes, whose value is given 
in Eq. (B.ll) of JUj- Such a canonical-momentum pseudofermion is related in [3j to the 
local av pseudofermion by a suitable Fourier transformation. The latter object occupies 
the sites of the effective av lattice [HUE]. Except for the discrete momentum values, 
the above pseudoparticle and pseudofermion have the same properties. Thus, all the 
energy eigenstates that span the PS can be described by occupancy configurations of 
av pseudofermions, —1/2 Yang holons, and —1/2 HL spinons jHIHj- The functional, 

Qtu(Qj) = 2n ®au, a >v>(qj,qf) AAWfeO , (2) 

ct'u',j 
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of equation ((TJ) was introduced in j3] and is such that Q® u (qj)/2 is found below to 
be an overall scattering phase shift. Here AN au (qj) = N au (qj) — N® u (qj) is the au 
branch bare-momentum distribution-function deviation relative to the ground state 
value and 7r $ ai/) a i v i (q, q') is defined in [3] and is found below to be an elementary two- 
pseudofermion phase shift. Note that Q%,{Qj) = for the initial ground state and thus 
cjj = qj for that state. 

Each transition from the initial ground state to a PS excited energy eigenstate can 
be divided into two elementary processes. The first process is a scattering-less finite- 
energy and finite-momentum excitation which transforms the ground state onto a well 
defined virtual state. This excitation involves the pseudofermion creation, annihilation, 
and particle-hole processes associated with the PS excited state and the discrete bare- 
momentum shift Qa U /L, whose possible values are 0, ±tt/L jHj, for au branches with 
finite occupancy in that state. For u > branches that excitation can involve a change 
in the number of discrete bare-momentum values. Although the au ^ cO, si branches 
have no finite pseudofermion occupancy in the initial ground state, one can define the 
values N* u = N^ u for the corresponding empty bands [TUl El H] • In this first step the 
pseudofermions acquire the excitation momentum and energy needed for the second- 
step scattering events. Thus, the virtual state is the in asymptote of the pseudofermion 
scattering theory. The second elementary step of the ground-state transition involves a 
set of elementary scattering events where all au pseudofermions or holes of momentum 
qj + Q° au / L of the in asymptote are the scatterers. Each of these elementary scattering 
events leads to a phase factor in the wave function of the au pseudofermions or holes 
given by, 

S a v,a>ui(q 3 , q 3 >) = e i27r fe>^') . (3) 

The scattering centers are the a'u' pseudofermions or holes of momentum q 3 > + Q au /L 
created in the ground-state - virtual-state transition and thus such that AN a i v i{qji) ^ 
0. Indeed, note that S aVj a'u'(qj,qj') = 1 for AN a > u i(qf) = 0. There is a one- 
to-one correspondence between the local rotated-electron occupancy configurations 
that describe the PS energy eigenstates and the local au pseudofermion occupancy 
configurations and —1/2 Yang holon and —1/2 HL spinon occupancies that describe 
the same states jlj. The corresponding effective au lattices have the same length 
L as the original lattice. Our analysis refers to periodic boundary conditions and 
the thermodynamic limit L — > oo. Under a ground-state - excited-energy-eigenstate 
transition, by moving the au pseudofermion or hole of initial ground-state momentum 
qj once around the length L lattice ring, its wave function acquires the following overall 
phase factor, 



N* 



S au (qj) = e tQ ° a » Yl II s au, a 'u'(qj,qf) 

a'u' j'=X 

- p iQuv(qj) . a — i o at* 

O , J 3 3 '"3 CXU ' 



(4) 
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Interestingly, Q a v(qj)/L is the net av pseudofermion or hole discrete canonical- 
momentum shift that arises due to the above transition [3J H] and thus in this equation, 



is such that Q a v(qj)/2 is a av pseudofermion or hole overall phase shift. Indeed, if when 
moving around the lattice ring the av pseudofermion or hole departures from the point 
x = —L/2 and arrives to x = L/2, one finds that linx^^ qx = qx + Qav(q)/2 
where q refers to the initial ground state. From Eqs. (0) and (J5J) it then follows 
that 7r & au ,a'u'(qj,qj') is an elementary two-pseudofermion phase shift. (If instead one 
considers x = and x = L, the overall phase shift and the two-pseudofermion phase 
shifts read Q au (q) and 2tt § au , a 'u'(qj, qj'), respectively [T3|. However, the choice of either 
definition is a matter of taste and the uniquely defined quantity is the above S matrix.) 

Several properties play an important role in the pseudofermion scattering theory. 
First, the elementary scattering processes associated with the phase factor (JHl) conserve 
the total energy and total momentum. Second, the elementary scattering processes are of 
forward-scattering type and thus conserve the individual in asymptote av pseudofermion 
or hole momentum and energy. These processes also conserve the av branch, usually 
called channel in the scattering language. Moreover, the scattering amplitude does 
not connect objects with different rj spin or spin. Last but not least, for each av 
pseudofermion or hole of initial ground-state momentum qj, the 5* matrix associated 
with the ground-state - excited-energy-eigenstate transition is simply the phase factor 
given in Eq. (jlj). For each excited energy eigenstate (out asymptote) the number of av 
pseudofermions plus the number of av pseudofermion holes whose S matrix is of the 
form flU) is given by N a + + Ea^cQ )S i 0(|AiV ai/ |) N* v . Here 6{x) = 1 for x > and 
9(x) = for x = 0. 

Importantly, the form of the scattering part of the overall phase shift (jnj), Eq. (J2J), 
reveals that the value of such a phase-shift functional is independent of the changes 
in the occupation numbers of the ±1/2 Yang holons and ±1/2 HL spinons. Thus, 
these objects are not scattering centers. Moreover, they are not scatterers, once their 
momentum values remain unchanged under the ground-state - excited-energy-eigenstate 
transitions. In turn, the pseudofermions and holes are scatterers and scattering centers. 
Since the cO pseudofermion is a ?7-spin-less and spin-less object and for v > the av 
pseudofermions are //-spin (a = c) and spin (a = s) singlet 2z/-holon and 2z/-spinon 
composite objects, respectively, their S matrix has dimension one: it is the phase factor 
(J3J). The factorization of the Bethe-ansatz (BA) bare S matrix for the original spin 
1/2 electrons is associated with the so called Yang-Baxter Equation (YBE) [Zj. On 
the other hand, the factorization of the S matrix (@J) in terms of the elementary S 
matrices S au , a > u >(qj,qji), Eq. (jH}, is commutative. Such a commutativity is stronger 
than the symmetry associated with the YBE and results from the elementary 5* matrices 
S a u,a'u'(qj,qj') being simple phase factors, instead of matrices of dimension larger than 
one. This seems to be inconsistent with all PS energy eigenstates being described 
by occupancy configurations which, besides cO pseudofermions, involve finite spin 1/2 
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spinons and 77-spin 1/2 holons ^U]- Indeed, the S matrix of finite 77-spin or spin 
objects has dimension larger than one. However, due to the correlations the quantum 
liquid self organizes in such a way that the scatterers and scattering centers are the cO 
pseudofermions, 77-spin singlet 2i/-holon composite cv pseudofermions, and spin singlet 
2i/-spinon composite sv pseudofermions. 

Let us clarify how the av pseudofermion S matrix (|1J) controls the unusual spectral 
properties of the model. Consider a av pseudofermion of canonical momentum q and a 
a'v' pseudofermion of canonical momentum q' such that the values q and q' correspond 
to a PS excited energy eigenstate and the initial ground state, respectively, and thus 
q' = q'. Importantly, from the use of Eq. it is found that the pseudofermion 
anticommutation relations introduced in [3] can be expressed solely in terms of the 
difference [q — q'] and the S matrix of the av pseudofermion associated with the excited 
state, 



{ftavi fq',a'u'} ~ $a 



q,a.V J 0£i ocV } w au,a'W 

1/2 



1/2 , ,„ Ira 

x 



N. 



01' 



sm([q - q'\/2) 

and the anticommutators between two creation or annihilation operators vanish. This 
reveals that the S matrix (0J) fully controls the pseudofermion anticommutators. Since 
within the PDT these anticommutators determine the value of the matrix elements 
between energy eigenstates jlj, it follows that the S matrix (jjj) controls the spectral 
properties. If it had dimension larger than one, the pseudofermion algebra would be 
much more involved, for the pseudofermion anticommutators would also be matrices of 
dimension larger than one. 

In reference jS] the excited states generated from the n = 1 and ra = ground 
state were described in terms of ±1/2 holon and ±1/2 spinon occupancy configurations. 
Following the analysis of Refs. [TT| IT2] for the related spin 1/2 isotropic Heinsenberg 
chain, the holes of the BA length-one spin string spectrum (spin singlet two-spinon 
composite si pseudoparticle spectrum) were identified in jS] with the spinons. Inspired 
in such an interpretation, the studies of the latter reference identified the holons with 
the holes of the BA distribution of k's |2 HO] spectrum [cO pseudoparticle spectrum]. 
This is behind the charge ±e found for the ±1/2 holons in |Hj, which is half of the value 
found in JUJ. However, the cO pseudoparticle and hole band occupancy configurations 
do not correspond to 77-spin SU (2) irreducible representations. Indeed, in (THj it is shown 
that for the whole Hilbert space all such representations exactly correspond to the BA 
charge string and ±1/2 Yang holon occupancy configurations. Following directly the 
analysis of Refs. ^], the studies of [H] consider that the ±1/2 holons and ±1/2 
spinons are the scatterers and scattering centers. This leads to two 4x4 5 matrices for 
holons and spinons, respectively, and a related 16 X 16 S matrix for the full scattering 
problem. In spite of being mathematically elegant and obeying the YBE, these matrices 
are not suitable for the description of the spectral properties. Moreover, provided that 
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within the x = and x = L boundary conditions one defines the overall phase shift 
as Qaviq), the phase shifts given in Eqs. (5.19)-(5.21) of the phase shifts given in Eqs. 
(5.19)-(5.21) of the first paper of |Hj, which appear in the entries of these matrices, are 
nothing but very particular cases of 77-spin-less and spin-less cO pseudofermion hole or 
spin-zero si pseudofermion hole overall phase shifts given in Eq. (0) JSj- Indeed, these 
phase shifts correspond to the n — 1 and m = initial ground state and the specific 
excited states considered in Refs. [7J |S]. Let qi or q\ be the bare momenta of the 
scattered cO or si pseudofermion hole, respectively, of the latter states. For the ?7-spin 
triplet, ?7-spin singlet, and 77-spin and spin doublet excited states considered in these 
references, it is found that n + Qco(qi) equals the phase shift 5ct and 5c s given in Eq. 
(5.19) and S v s in Eq. (5.21) of the above paper, respectively. For the spin triplet, spin 
singlet, and spin and ^-spin doublet excited states, Q s i(q'i) equals the shift functions 
5st and Sss given in Eq. (5.20) and Ss v given in Eq. (5.21), respectively [12]. Thus, the 
BA phase shifts of Refs. [HE] are particular cases of the cO and si pseudofermion hole 
overall phase-shift functionals of Eq. (jHJ) and are associated with a set of excited states 
which span a subspace smaller than the PS of the one- and two-electron excitations. 
According to the studies of (THj, for all the transitions associated with these phase shifts 
the deviations in the 77-spin and spin values are provided by the ±1/2 Yang holons and 
±1/2 HL spinons occupancy changes, respectively, which do not contribute to the phase- 
shift values. In turn, the holes created in the cO and si bands by these transitions are 
both scatterers and scattering centers and it follows from the analysis of ^U] that they 
do not correspond to single ±1/2 holons and ±1/2 spinons, respectively. Moreover, the 
phase shifts of Refs. jZHH] were evaluated up to an overall constant term by the method 
of [Hj. Equation © provides the full phase shift value and reveals that the above extra 
7r in Sqt, 8cs, and S v s is not physical, as discussed elsewhere [T3]. The results reported 
here also apply to other models. For instance, for the isotropic Heisenberg chain it is 
found that the phase changes for the spin singlet and triplet excited states given in Eq. 
(11) of (TTj equal the phase shifts for the same states of a scattered hole of the zero-spin 
two-spinon si pseudofermion spectrum. Thus, for the study of the spectral properties, 
these two excited states correspond to two si pseudofermion hole S matrices, rather 
than to the single 4 x 4 S matrix of Eq. (5.1) of ^2j. This analysis can also be extended 
to the same model with an odd number of lattice sites. 

While, through the anticommutators the use of the S'-matrix introduced here 
leads to a successful description of the spectral features observed in real materials [Sj, 
the 16 x 16 S matrix of [H] is unsuitable for such a task. Indeed, independent ±1/2 
holons and ±1/2 spinons that are not part of composite pseudofermions are neither 
scatterers nor scattering centers. Interestingly, these objects remain invariant under 
the electron - rotated-electron unitary, whereas the pseudofermion and holes are not in 
general invariant under such a transformation. 

The method for evaluation of the finite-energy spectral-weight distributions of a 
ID correlated metal introduced in [3] fully relies on the scattering theory introduced 
here. The exotic metallic quantum phase of matter found for quasi-lD compounds [3] 
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by use of such a method is expected to emerge at finite energies in carbon nanotubes, 
ballistic wires, and systems of cold fermionic atoms in one- dimensional optical lattices 
with on-site atomic repulsion [Oj. This confirms the general scientific interest of the 
scattering theory introduced here. While in this Letter it is applied specifically to 
the ID Hubbard model, the theory is of general nature for many integrable quantum 
problems and therefore will have wide applicability. 
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